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ABSTRACT 



O ' The spectral energy distributions (SEDs) of dust disks are widely used to infer 

dust properties (compositions and sizes) and disk structures (dust spatial distributions) 
which might be indicative of the presence or absence of planets or smaller bodies like 
, asteroids and comets in the disk. Based on modelling of the SED of HR4796A, a young 

main-sequence star with the largest fractional infrared (IR) emission, we show that the 
SED alone is not a good discriminator of dust size, spatial distribution (and composition 
if no spectroscopic data are available). A combination of SED, mid-IR spectroscopy, 
and coronagraphic near-IR imaging of scattered starlight and mid-IR imaging of dust 
■ thermal emission provides a better understanding of these properties. 

Subject headings: circumstellar matter — dust, extinction — infrared: stars — plane- 
tary systems: protoplanetary disks — stars: individual (HR 4796A) 



1. Introduction 

Over the past 2 decades, impressive evidence has been assembled for the existence of circumstel- 
lar dust disks around main-sequence (MS) stars as well as pre-MS stars (T Tauri stars and Herbig 
Ae/Be stars), post-MS stars (red giants), and a white dwarf (see Zuckerman 2001 for a review). 
A wide variety of observational techniques have been employed to study the formation/evolution 
and physical/chemical properties of dust disks: optical and near- infrared (IR) imaging of scattered 
stellar light, photometric measurements of dust thermal emission from near-IR to submillimeter, 
spectroscopic observations of mid-IR dust emission features and gas emission lines and ultraviolet 
(UV) and visible gas absorption lines. The spectral energy distribution (SED) is of particular 
interest in inferring the size and composition of dust grains and the disk structure (dust spatial 
distribution). However, the limitations of this method have not been adequately explored; for 
example, in modelling the (3 Pictoris SED, Li & Greenberg (1998) found that the dust spatial 
distribution is coupled with the distribution of dust sizes, i.e., the distribution of dust in the disk 
and dust sizes cannot be uniquely determined simultaneously by the SED alone. 
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It is the aim of this Letter to quantify the limitations on the information derived from the 
SED modelling concerning the dust properties (composition and size) and disk structure. We 
take the SED of HR4796A, a nearby (distance to Earth d ~ 67 ± 3pc) young MS star (age 
8±3Myr) of spectral type AO V (effective temperature T e g- « 9500 K) for comparison with model 
results. HR4796A has the largest fractional IR luminosity relative to the total stellar luminosity 
(Lir/L* «5x 10 -3 ) among the ~ 1500 A-type MS stars in the Yale Bright Star Catalogue (Jura 
1991). The dust disk in orbit around HR4796A has been extensively studied, both observationally 
and theoretically (see Zuckerman 2001 and references therein). With its relatively well determined 
dust and disk properties, HR4796A serves as a good comparison basis for the SED modelling efforts 
described in this Letter. We stress that the main purpose of this Letter is not to carry out a detailed 
study of the HR4796A dust disk which we defer to a subsequent paper (A. Li & J.I. Lunine 2002, 
in preparation). We will first outline our approach in §2. We then discuss in §3 the degeneracy 
between the dust spatial distribution and dust sizes under the assumption of pure silicate dust. In 
§4 we show that the dust composition is not well constrained by the observed SED unless mid-IR 
spectroscopic dust emission features are available. We discuss in §5 the possible dust composition 
and morphology from the evolutionary point of view that circumstellar dust disks around (pre- 
)MS stars are formed through the coagulation of interstellar solids. We also summarize our major 
conclusions in S5. 



2. Modelling the Dust IR Emission 

Grains in the optically thin dust disk of HR4796A absorb the stellar UV/optical photons and 
then re-radiate them in the IR. Very small grains or large molecules (with spherical radius smaller 
than 100 A) which are subject to single-photon heating (Draine & Li 2001) will not be considered 
here although their presence in the disk can not be ruled out (e.g., polycyclic aromatic hydrocarbon 
molecules are seen in the disk of HD 100546 through their 3.3, 6.2, 7.7, 8.6, and 11.3 fim emission 
features [Malfait et al. 1998]). For grains > 100 A we first use Mie theory to calculate absorption 
cross sections assuming a spherical shape. We then calculate grain equilibrium temperatures (and 
corresponding IR emission) by balancing absorption and emission. For a given dust size distribution 
and a given disk structure (dust spatial density distribution), the emergent IR emission spectrum 
can be obtained by integrating over the dust size range, and over the entire disk. 

We have the following parameters to be specified or constrained: (1) dust composition - we will 
consider amorphous silicate, amorphous carbon, and cometary grains (porous aggregates of small 
silicate and carbon dust; Greenberg 1982); (2) dust sizes - we will consider a power law dust size 
distribution dn(a)/da oc a~ a characterized by the lower-cutoff a m i n , upper-cutoff a max and power- 
law index a (where a is the spherical radius); (3) dust spatial density distribution - we will consider 
two dramatically different distribution functions: a power law distribution dn(r)/dr oc r~@ and a 
Gaussian distribution dn{r)/dr oc exp[— 41n2{(r — ro)/A} 2 ], the former is characterized by the 
disk's inner boundary r; n , outer boundary r out , and power-law index (3; the latter is characterized 
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by the radial position ro where dn(r)/dr peaks and the full width half maximum (FWHM) A. 

Juraet al. (1995, 1998) and Augereau et al. (1999) estimated that grains with radius exceeding a 
few micron are stable against radiation pressure. Telesco et al. (2000) estimated the "characteristic" 
radius for the 10-20 /jm mid-IR-emitting grains to be ps 1 — 1.5 /im. We therefore adopt a m \ n = 1 [xm. 
We take a max = 1 cm (this is a noncritical parameter since grains larger than ~ 100 are like 
blackbodies and their IR emission spectra are size-insensitive). 

To be conservative, we first set the inner boundary at r in = 0.15 AU inside of which micron- 
sized silicate and carbonaceous grains would be heated to > 1500 K and evaporate. The outer 
boundary is taken to be r out = 250 AU which is expected from the disk truncation caused by the 
tidal effects of HR4796B, a companion star of HR4796A (Jayawardhana et al. 1998). Other values 
for rj n and r out estimated from the near- and mid-IR imaging observations will also be discussed 
(see §3). For the Gaussian-type dust spatial distribution, we take the peak distance (from the 
central star) of the distribution to be ro = 70 AU as reflected from the near-IR imaging of scattered 
starlight (Schneider et al. 1999) and mid-IR imaging of dust thermal emission (Jayawardhana et al. 
1998; Koerner et al. 1998; Telesco et al. 2000). Therefore, we are only left with two free parameters: 
the dust size distribution power index a and the dust spatial distribution power index (3 or the 
FWHM A of the Gaussian-type spatial distribution. 



3. Dust Spatial Distributions and Dust Sizes 

We first consider compact silicate grains. We adopt the indices of refraction of the Draine & 
Lee (1984) "astronomical" silicate. We approximate the HR4796A radiation field by the Kurucz 
model atmosphere spectrum for AO V stars (Kurucz 1979). We will compare our model results with 
the available photometric data compiled by Augereau et al. (1999) for the HR4796A disk. 

Assuming a power-law dust spatial distribution for r- ul = 0.15 AU and r out = 250 AU, the best 
fit to the observed SED is provided by a w 3.8 and (3 ps —1.6 (see Figure 1) - the density of 
dust increases outward from the star through the whole disk. Since various studies suggest that 
the HR4796A disk has an inner hole at r ~ 40 — 60 AU and an outer edge sharply truncated 
at ~ 80 — 130 AU (Jura et al. 1993; Jura et al. 1995; Jayawardhana et al. 1998; Koerner et al. 
1998; Schneider et al. 1999; Wyatt et al. 1999; Telesco et al. 2000), we also consider models with 
r; n = 40 AU and r ou t = 130 AU. To obtain a satisfactory fit to the observed SED, an even steeper 
outward increase of dust distribution is needed (for enhancing the amount of warm dust). In Figure 
1 we also present the best-fit model spectrum for r- m = 40 AU and r out = 130 AU obtained from 
a ss 3.6 and [3 pa —2.7. We have also modelled the HR4796A SED in terms of broken power- 
laws for both the dust spatial distribution and dust sizes: dn(r)/dr oc r _/3l ; dn(a)/da oc a~ ai for 
r in ^5 f ^ tq and dn(r)/dr oc r~@ 2 , dn(a)/da oc a~ a ' 2 for ro < r < r out . These 2-power-law models 
also require j3 < and we do not see significant improvements compared with single-power-law 
models. 
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Motivated by the NICMOS (Near-IR Camera and Multi-Object Spectrometer) discovery that 
the surface brightness of the HR4796A disk's scattered light sharply peaks at ~ 70 AU (Schneider 
et al. 1999), we adopt a Gaussian function peaking at ro = 70 AU for the dust spatial distribution. 
As long as the FWHM A < 30 AU, the bulk of the dust lies at 40 < r < 130 AU and we expect 
little difference between the [r in = 0.15AU,r out = 250 AU] model and the [r^ n = 40AU,r out = 
130 AU] model. Therefore, for the Gaussian-type distribution, we only consider r- m = 0.15 AU and 
r out = 250 AU. Figure 1 plots the best-fit spectrum calculated from r- m = 0.15 AU, r out = 250 AU, 
ro = 70 AU, A = 20 AU, and a = 3.5. This appears to be consistent with the NICMOS detection 
of a sharply bounded and narrow (A < 17 AU) ring-like disk (Schneider et al. 1999). 

It is seen in Figure 1 that, although the dust spatial distributions are strikingly different, 
all three models are able to provide similarly good fits to the observed SED except the IRAS 
[Infrared Astronomical Satellite] 60 /im data (we note that the IRAS photometric uncertainty given 
by Augereau et al. (1999) might have been underestimated [Beichman et al. 1989]), provided that 
the dust sizes are different. 

4. Dust Composition 

We have seen in §3 that pure solid silicate grains with various sizes and spatial distributions are 
successful in reproducing the SED of HR4796A. Similarly, the HR4796A SED can also be fitted by 
pure solid amorphous carbon grains using the index of refraction of Rouleau & Martin (1991). For 
illustration, we present in Figure 1 the best-fit single-power-law model spectrum calculated from 
ri n = 0.15 AU, r out = 250 AU, f3 = —5.0, and a = 4.2. Although it is unlikely that (proto)planetary 
disk dust is mainly made of carbonaceous material, the observed SED alone is unable to rule out 
the pure carbon dust model unless the 9.7/xm Si-0 and 18 /im O-Si-0 features are detected. 1 We 
will discuss this further in §5. 

Jura et al. (1998) argued that the HR4796A disk are composed of cometary icy grains. The 
HR4796A SED was closely reproduced by Augereau et al. (1999) in terms of a cold interstellar 
dust component and a hot cometary dust component. In this Section we also model the HR4796A 
SED by cometary grains. Following Greenberg (1982), we model cometary dust as fluffy aggregates 
of interstellar silicate and carbonaceous grains. 2 We take the grain porosity (the volume fraction 



lr The 8-13 /im spectroscopic observations of the HR4796A disk show that the thermal emission in the silicate 
feature is very weak (Sitko, Lynch, & Russell 2000). 

2 In the framework of the Greenberg comet model (Greenberg 1982, 1998; Greenberg & Li 1999), the silicate 
dust and carbon dust are physically related through a core-mantle structure. In this work the physical relationship 
between the silicate dust and carbon dust does not matter since the Bruggeman effective medium theory employed to 
calculate dust optical properties does not distinguish inclusions from matrix (Bohren & Huffman 1983). Therefore, 
the Greenberg comet model is not just limited to the core-mantle interstellar dust model (Li & Greenberg 1997), 
but also applicable to other popular dust models such as the separate silicate/graphite model (Mathis, Rumpl, & 
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Fig. 1. — Comparison of the observed spectral energy distribution of the HR4796A dust disk to theoretical IR 
emission spectra calculated from (1) heavy solid line - pure compact silicate grains (mass density p B u ~ 3.5 gem -3 ) 
with a power law size distribution dn(a)/da oc a~ a where a — 3.8 and a G [1 /im, 1cm], and a power-law spatial 
distribution dn(r)/dr oc r"* 3 where /3 = -1.6 (r G [0.15 AU, 250 AU]) and total dust mass m dust « 4.98 x 10 26 g; 
the slope change at A ~ 7 /im is due to the onset of the 10 /im Si-O silicate vibrational mode; the little bump at 
10 fim is the Si-O silicate band; (2) dotted line - silicate grains with a flatter power-law size distribution of a — 3.6 
(a G [1 /im, lem]) and a steeper and less extended spatial distribution of /3 = —2.7 (r G [40 AU, 130 AU]) and 
Wdust « 4.53 x 10 26 g; (3) dashed line - silicate grains with a power-law size distribution of a = 3.5 (a G [1 /im, 1 cm]) 
and a Gaussian spatial distribution of dn(r)/dr oc exp[— 41n2{(r — ro)/A} 2 ] with its peak at ro = 70 AU and a 
FWHM A = 20 AU (r 6 [0.15 AU, 250 AU]) and m dus t ~ 3.46 x 10 26 g; (4) thin solid line - pure compact amorphous 
carbon grains (mass density p am c ~ 1.8 gem -3 ) with a power law size distribution of a = 4.2 (a £ [1 /mi, 1cm]) and 
a power-law spatial distribution of /3 = -5.0 (r G [0.15 AU, 250 AU]) and m dust w 1.89 x 10 26 g. 

of vacuum) to be 0.6 (Augereau et al. 1999). The volume ratio of the silicate component to the 
carbonaceous component is taken to be 1.0 as derived from the in situ measurement of the coma 
dust of comet Halley (Kissel & Krueger 1987). We use Mie theory and the Bruggeman effective 
medium theory (Bohren & Huffman 1983) to calculate the absorption cross sections for porous 
cometary grains. Similar to the pure solid silicate dust model (see §3), the cometary dust model 
is also able to provide reasonably good fits to the HR4796A SED with various dust sizes and dust 
spatial distributions. In Figure 2 we show the model spectra calculated from cometary grains with a 
power-law size distribution and a power-law or Gaussian spatial distribution: (1) a = 3.9, [3 = —2.0 



Nordsieck 1997, Draine & Lee 1984, Weingartner & Draine 2001, Li & Draine 2001) and the composite dust model 
(Mathis 1996). However, in the dense protostellar environment, it is possible that graphite grains may be destroyed 
by chemical attacks of O, H2O (Draine 1985). 
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and re [0.15 AU, 250 AU]; (2) a = 3.6, = -3.7 and r € [40AU, 130AU]; (3) a = 3.5, r = 70 AU, 
A = 15 AU and r € [0.15 AU, 250 AU]. 
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Fig. 2. — Comparison of the observed SED of the HR4796A dust disk to theoretical IR emission spectra of cometary 
grains (composed of silicate, carbonaceous dust and vacuum): (1) heavy solid line - for a power law size distribution 
of a = 3.9 (a £ [1 /im, 1 cm]) and a power-law spatial distribution of (3 — —2.0 (r 6 [0.15 AU, 250 AU]) and total dust 
mass mdust ~ 3.23 x 10 26 g; (2) dotted line - for a power law size distribution of a = 3.6 (a g [1 /im, 1 cm]) and a power- 
law spatial distribution of /3 = —3.7 (r £ [40 AU, 130 AU]) and md us t ~ 4.47 x 10 26 g; (3) dashed line - for a power 
law size distribution of a = 3.5 (a £ [1 /im, 1 cm]) and a Gaussian spatial distribution of ro = 70 AU and A = 15 AU 
(r 6 [0.15 AU, 250 AU]) and m dus t ~ 3.00 x 10 26 g. In all models the volume fractions of the silicate, carbonaceous, 
and vacuum components are taken to be 0.2, 02., and 0.6, respectively (effective mass density fa 1.06 gem -3 ). Also 
plotted is the model spectrum (thin solid line) calculated from a Gaussian distribution (r- m ~ 0.15 AU, r ou t = 250 AU, 
ro = 70 AU, A = 15 AU) of grains with compositions derived in the context of the formation and evolution of dust 
disks (m C arb/Ti. s ii = 0.7 and a porosity of P = 0.6 for r < 30AU; compact icy grains with m car b/m s ii = 0.7 and 
micc/[m s ii + m carb ] = 0.8 for r > 30 AU; see §5) with a power-law size distribution of a = 3.6 and total dust mass 
m duat » 1.94 x 10 26 g. 



5. Discussions 

It is shown in §3 and §4 that models with various compositions, sizes, and spatial distributions 
are able to reproduce the observed SED of the HR4796A disk reasonably well. It was also shown 
by Jura et al. (1998) that icy grains with a typical radius near 100 /im are able to explain the 
HR4796A SED. Dust thermal emission depends on its absorption and emission properties which 
are determined by its size and optical properties. It would not be surprising for a wide range 
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of dust materials with properly chosen sizes to be able to fit the SED. However, we should not 
be too pessimistic: such spectrum will be useful when combined with other constraints on the 
composition of the dust in circumstellar disks around (pre-)MS stars. The coagulation of interstellar 
grains that results in fluffy and inhomogeneous aggregates occurs in cold, dense molecular clouds 
and protostellar and protoplanetary dust disks and plays an important role in the formation of 
planetary systems (Weidenschilling & Cuzzi 1993). We can therefore approximately derive the 
proportional composition of the dust in circumstellar disks from the abundances of the condensable 
elements (C, N, O, Si, Fe, and Mg), 3 assuming protostellar activities impose little modification on 
protostellar grain compositions (see Beckwith, Henning, & Nakagawa 2000). 

Let [X/H] be the cosmic abundance of X relative to H (we assume the cosmic elemental abun- 
dances are those of the solar values: [C/H] Q ~ 391 parts per million (ppm), [N/H] Q ~ 85.2 ppm, 
[O/H] » 545 ppm, [Mg/H] » 34.5ppm, [Fe/H] Q » 34.4ppm, and [Si/H] » 28.1ppm [Sofia & 
Meyer 2001]); [X/H] gas be the amount of X in gas phase ([C/H] gas « 140ppm, [N/H] gas « 61 ppm, 
[0/H] gas rs 310 ppm; Fe, Mg and Si are highly depleted in dust; see Li & Greenberg 1997 and 
references therein); [X/H] dust be the amount of X relative H locked up in dust ([C/H] dust = 
[C/H] -[C/H] gas « 251 ppm, [N/H] dust « 24.2 ppm, [0/H] dust « 235 ppm, [Mg/H] dust « 34.5 ppm, 
[Fe/H] dust ~ 34.4 ppm, [Si/H] dust ~ 28.1 ppm). Assuming a stoichiometric composition of MgFeSi04 
for interstellar silicates, the total mass of silicate dust per H atom is m si \ « [Fe/H] dust + 
[Mg/H] dust nug + [Si/H] dust ^si + [0/H] sil /i ~ 5.61 x 10~ 3 /i H where /U X is the atomic weight of X 
in unit of m « 1.66 x lO" 24 g, and [0/H] sil « 4 ([Fe/H] dust + [Mg/H] dust + [Si/H] dust )/3 « 129ppm 
is the amount of O in silicate dust per H atom (i.e., we assign 4 O atoms for the average of the 
Fe, Mg, and Si abundances). The carbonaceous dust component is dominated by C, with lit- 
tle H, N, and O (we assume H/C=0.5, O/C=0.1). The total mass of carbon dust per H atom 
is m carb « [C/H] dust fi c + [N/H] dust m + 0.5 [C/H] dust m + 0.1 [C/H] dust Mo « 3.88 x 10" 3 m . 
The C, O, and N atoms left over after accounting for the silicate and carbon dust components 
are assumed to condense in icy grains in the form of H2O, NH3, CO, CO2, CH3OH and CH4 
(following Greenberg [1998], we assume CO:CO 2 :CH3OH:CH4:H 2 CO=10:4:3:l:l). The total mass 
of icy grains per H atom is m lcc ~ m? e + mf ce + m™ tcr , where the mass of C-containing ice 
mp cc » [C/H] gas ^ c + [C/H] gas (22/xo + 18/x H )/19 » 2.87 x 10~ 3 m ; the mass of NH 3 ice mf ce » 
[N/H] gas ( MN + 3a*h) ~ 4.10 x 10-Vh; the mass of water ice m™ ter » [0/H] watcr (/x + 2^ H ) ~ 
4.12 x 10- 3 ^ H ; [0/H] watcr « [O/H] - [0/H] sil - 0.1 [C/H] dust - 22 [C/H] gas /19 « 229 ppm is the 
amount of O locked up in H2O ice (we assume H2O contains all the remaining available O). 

Therefore, as a first approximation, we may assume a mixing ratio of m caiT ^/m si \ ps 0.7 and 
"ike/ (m s ii + m car b) ~ 0.8 for cold regions (for hot regions where ices sublimate the dust can be 
simply modelled as porous aggregates of silicate and carbon particles with m car b /m s n ~ 0.7). 
This does not deviate much from the in situ measurements of cometary dust (m car b / 'm s ii ~ 0.5, 
^icc/b^sii + "^carb] ~ 1-0; see Greenberg & Li 1999 and references therein) which is often suggested 



3 Some H will be present, mostly in combination with O, C, and N. 
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as porous aggregates of unaltered interstellar dust (Greenberg 1982; Greenberg & Li 1999). The 
porosity is a free parameter ranging from that of diffuse cloud interstellar dust (~ 0.45, Mathis 
1996) to that of very fluffy cometary dust (> 0.9, Greenberg & Li 1999). For illustration, we plot in 
Figure 2 the model spectrum calculated from a Gaussian distribution (n n = 0.15 AU, r out = 250 AU, 
r = 70 AU, A = 15 AU) of grains with (1) a power law size distribution (a min = 1 /xm, a max — 1 cm, 
a = 3.6), (2) mcarb/w-sii = 0.7 and a porosity of P = 0.6 for r < 30 AU; and (3) compact icy grains 
with m car b/m s ii = 0.7 and m,i cc /(m s ii + m car b) = 0.8 (porous grains of P = 0.6 become compact 
after filled with ices of an amount of mi cc /[m s i\ + m car b] = 0.8) for r > 30 AU. This will be discussed 
in detail in a subsequent paper (A. Li & J.I. Lunine 2002, in preparation). 

The fact that pure amorphous carbon grains are also able to account for the observed SED (see 
§4) reinforces the importance of combining observations with theoretical calculations of dust com- 
position in the context of the formation and evolution of dust disks. We note that the non-detection 
of the silicate emission features in the HR4796A disk (Sitko et al. 2000) does not necessarily imply 
the predominance of non-silicate dust in the disk. It may just imply the lack of small and hot 
silicate grains. On the other hand, the contradistinction between the various dust spatial distribu- 
tions, which all provide a reasonably good fit to the observed SED (see §3 and §4), indicates the 
importance of direct disk imaging. 

In summary, the spectral energy distributions of dust disks alone are not necessarily able 
to constrain the dust compositions, sizes, and spatial distributions. The dust spatial and size 
distributions are coupled. Caution should be taken in discussing the presence/absence of planets, 
comets, and asteroids in the disk solely based on the observed SED. We argue that grains in 
circumstellar disk around (pre-)MS stars are composed of silicate, carbonaceous dust (and ices in 
cold regions) and vacuum with a mixing ratio of m car b /m s \\ f» 0.7 and mi cc /(m si i + m car b) ~ 0.8. 
A combination of compositional considerations, SED, mid-IR spectroscopy, coronagraphic near-IR 
imaging of scattered starlight and mid-IR imaging of dust thermal emission will allow us to better 
understand the properties of circumstellar dust disks. 

We dedicate this paper to Prof. J. Mayo Greenberg who passed away on November 29, 2001. A. 
Li thanks the University of Arizona for the "Arizona Prize Postdoctoral Fellowship in Theoretical 
Astrophysics". This research was supported in part by NASA grant NAG5-10450. 
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